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Abstract 

Cooperative communication has been shown to provide significant increase of transmission re¬ 
liability and network capacity while expanding coverage in cellular networks. The present work is 
devoted to the investigation of the end-to-end performance and power allocation of a maximum- 
ratio-combining based regenerative multi-relay cooperative network over non-homogeneous scattering 
environment, which is the case in realistic wireless communication scenarios. Novel analytic expressions 
are derived for the end-to-end symbol-error-rate of both M—ary Phase-Shift Keying and M—ary 
Quadrature Amplitude Modulation over independent and non-identically distributed generalized fading 
channels. The offered results are expressed in closed-form involving the Lauricella function and can 
be readily evaluated with the aid of a proposed computational algorithm. Simple expressions are also 
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derived for the corresponding symbol-error-rate at asymptotically high signal-to-noise ratios. The derived 
expressions are corroborated with respective results from computer simulations and are subsequently 
employed in formulating a power optimization problem that enhances the system performance under 
total power constraints within the multi-relay cooperative system. Furthermore, it is shown that optimum 
power allocation provides substantial performance gains over equal power allocation, particularly, when 
the source-relay and relay-destination paths are highly unbalanced. 

Index Terms 

Asymptotic analysis, decode-and-forward, digital modulations, generalized fading channels, maximum- 
ratio combining, multi-relay systems, optimum power allocation. 

I. Introduction 

Cooperative transmission methods have attraeted signifieant interest over the past deeade due 
to their applicability in size, power, hardware and price constrained devices such as cellular 
mobile devices, wireless sensors, ad-hoc networks and military communications [Hl-llll. Such 
systems exploit the broadcast nature and the inherent spatial diversity of wireless paths and are 
typically distinguished between regenerative (decode-and-forward) or non-regenerative (amplify- 
and-forward) relaying schemes. In general, the digital processing nature of regenerative relaying 
is considered more efficient than non-regenerative relaying, since the latter typically requires 
costly RF transceivers in order to avoid forwarding a noisy version of the signal ifTOll - lfT^ . 

The performance of cooperative systems can be substantially improved by optimum allocation 
of the limited overall power to the source and relays of the network in order to minimize the 
overall energy consumption for given end-to-end performance specifications. Among others, this 
can be efficiently achieved by accurately accounting for the detrimental effects of multipath 
fading lfT7]l - ll38]l and the references therein. Based on this, the authors in ll^ derived upper 
and lower bounds for the outage probability (OP) of multi-relay decode-and-forward (DF) 
networks over independent but non-identically distributed (i.n.i.d) Nakagami—m fading channels. 
In the same context, the authors in [|4^ analyzed the symbol error rate (SER) and OP of DF 
systems with relay selection over i.n.i.d Nakagami—m fading channels, with integer values of 
m, whereas a comprehensive analytical framework for a dual-hop multi-antenna DF system 
under multipath fading was derived in iHTll . Likewise, the performance of DF systems over 
different fading environments was investigated in [I^ - [|43]| whereas analysis for the SER of 
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dual-hop DF relaying for M-ary phase-shift keying (M—PSK) and M-ary quadrature amplitude 
modulation (M—QAM) over Nakagami—m fading ehannels was reported in ll46ll . In addition, 
optimum power alloeation (OPA) in dual-hop regenerative relaying with respeet to pre-defined 
thresholds for the SER and OP was analyzed in Il4^ and Il471l . respeetively. This problem was 
also addressed in [|4^ and ll49l for multi-node DF relaying based on asymptotie SER and for 
a given network topology over Rayleigh fading ehannels, respeetively. Finally, the authors in 
[ISOll proposed power alloeation sehemes for the ease of multi-relay DF eommunieations in the 
high-SNR regime over Nakagami—m fading ehannels. 

Nevertheless, all reported investigations have been earried out in the eontext of either asymp¬ 
totie or dual-hop seenarios as well as eonsidering only Rayleigh or Nakagami—m fading ehan- 
nels, whieh is mainly due to the presenee of eumbersome integrals that involve eombina- 
tions of elementary and speeial funetions IISTII - IISTII and the referenees therein. However, it 
is reealled that these fading models are based on the underlying eoneept of homogeneous 
seattering environments, whieh is not praotieally realistie sinee surfaees in most radio propagation 
environments are spatially eorrelated lISSll . This issue was addressed in [|59]| by proposing the 
Tj—n distribution, whieh is a generalized fading model that has been shown to provide partieularly 
aeeurate fitting to realistie measurement results, while it ineludes as speeial eases the well known 
Rayleigh, Nakagami—m and Hoyt distributions Il59ll - ll63 . Based on this, several eontributions 
have been devoted to the analysis of various eommunieation seenarios over generalized fading 
ehannels that follow the t] — distribution, see, e.g., If62]| - ll68ll and the referenees therein. 
Motivated by this, the present work is devoted to the evaluation of the end-to-end SER in 
regenerative eooperative eommunieation systems with multiple relays for M—PSK and M—QAM 
eonstellations over generalized fading ehannels as well as the eorresponding OPA analysis. 
Speeifieally, the eontributions of this work are summarized below: 

• Exaet elosed-form expressions are derived for the end-to-end SER of M—PSK and M—QAM 
based multi-relay regenerative networks over generalized multipath fading environments for 
both i.n.i.d and i.i.d seenarios. 

• Simple asymptotie expressions are derived for the above seenarios for high SNR values. 

• The eorresponding amount of fading is derived for quantifying the respeetive fading severity. 

• Optimal power alloeation based on the eonvexity of the derived asymptotie expressions is 
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Fig. 1. Multi-node dual hop cooperative relay system. 


formulated, to minimize the eorresponding SER under sum-power eonstraint in all nodes. 

• The derived expressions are employed in evaluating the performanee of the eonsidered 
system and extraeting useful insights. 

• It is shown that a 3dB gain is aehieved by OPA for even a small number of relays. 

• It is shown that post-Rayleigh fading eonditions result to an improved performanee by up 
to 4dB eompared to eommunieations over severe fading ehannels. 

• It is shown that a maximum gain of about 21dB oeeurs, eompared to ordinary direet 
eommunieation, even if only few nodes are employed in non-severe fading eonditions. 
This renders the resouree eonstrained eommunieation system a meaningful alternative for 
inereasing the quality of serviee of demanding emerging wireless systems. 

• A simple algorithm is proposed for the eomputation of the generalized Laurieella funetion. 

The reminder of the present paper is organized as follows: Seetion II revisits the eonsidered 

system and ehannel models. The exaet and asymptotie SER expressions for M—QAM and 
M—PSK eonstellations over generalized multipath fading ehannels are derived in Seetion III. 
The optimum power alloeation seheme based on sum-power eonstraint is provided in Seetion 
IV while Seetion V presents the eorresponding numerieal results along with related diseussions 
and insights. Einally, elosing remarks are given in Seetion VI. 

II. System and Channel Models 

A. System Model 

We eonsider a multi-node eooperative radio aeeess system eonsisting of a souree node S, 
intermediate relay nodes Rk, with /c = {l,2,---,iT}, and a destination node D, as depieted in 
Eig.l. Eaeh node in the system is equipped with a single antenna while a half-duplex deeode- 
and-forward protoeol is adopted. Eurthermore, in order to avoid inter-relay interferenee the nodes 
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in the system transmit signals through orthogonal ehannels. In phase I, the souree broadeasts 
the signal to the destination and to all relay nodes in the network whieh yields 

ys,D = V^o:s,dX + ns^D (1) 


and 


ys,Rk 



( 2 ) 


where Pq is the souree transmit power, x is the transmitted symbol with normalized unit energy 
in the first transmission phase, as,D and as,R^. are the eomplex fading eoeffieients from the 
souree to the destination and from the souree to the relay, respeetively, whereas us^d and 
ris^R^ represent the eorresponding additive-white-Gaussian noise (AWGN) with zero mean and 
varianee Nq. In the next time slot, if the relay deeodes eorreetly, then it forwards the deeoded 
and re-eneoded signal to the destination with power Pr^ = Pr^; otherwise, Pr^ = 0 yielding 

yRk,D = -sj PRk(^Rk,DX + nR^,^D (3) 

where £> denotes the eomplex fading eoeffieient from the P'^ relay to the destination and 
riR^ R is the eorresponding AWGN. It is also assumed that eaeh path experienees narrowband 
multipath fading that follows the y — distribution and that MRC diversity seheme is employed 
at the destination. As a result, the eombined output reeeived signal is expressed as follows 

K 

yo = Woys,D + '^kyRt,,D (4) 

k=l 

where Wq = y/P^a*g j^/Nq and Wk = PR^a\^ d/^o denote the optimal MRC eoeffieients for 
ys,D and yR^^^o, respeetively with (•)* representing the eomplex eonjugate operator. 


B. Generalized Multipath Fading Channels 

It is reealled that y—p distribution has been shown to aeeount aeeurately for small-seale 
variations of the signal in non-line-of-sight eommunieation seenarios. This fading model is 
deseribed by the two named parameters, y and p, and it is valid for two different formats 
that eorrespond to two physieal models [[59ll . The probability density funetion (PDF) of the 
instantaneous SNR y = \ a \'^ P/Nq is given by 


A (7) = 


2y/^ /i^+2 2 


■ exp 


2pyh\ (2pHy\ 


7 / ^ V 7 / 


(5) 
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TABLE I 

Relation Between Distribution and Other Common Fading Distributions. 


Fading Distribution 

Format-1 

Format-2 

V - P 

h = (1 + 77-7 + 77)/4, H = (r;"7 

h = 1/(1 -7,2), H = 7,/(l - 7,2) 

Nakagami— 77T 

^ = m, r) —> 0 or r) 00 

77 —>■ ±1 

^ = ml2, r) 1 

77 0 

Nakagami —g (Floyt) 

/T = 0.5, r] = 

E = - '))/(! + v) 

Rayleigh 

fi — 0.5 ,77 — 1 

= 0.5 ,77 = 0 


where 7 = E{j) = {P/Nq)Q is the average SNR per symbol, with E{-) denoting statistieal 
expeetation, and = E{\ a p) denotes the channel variance. Furthermore, r(-) and /„(•) 
denote the Euler gamma function and the modified Bessel function of the first kind, respectively 
ll69l . The parameters h and H are given in terms of 7 in two formats, as depicted in Table I 
along with the fading distributions that are included in 77 —/i as special cases. In terms of physical 
interpretation 7 denotes the scattered-waves power ratio between the in-phase and quadrature 
components of the scattered waves in each multipath cluster in Format-1 and the correlation 
coefficient between the in-phase and quadrature components of each multipath cluster in Format- 
2. Fikewise, /i = £'^( 7 ) (1 -f {El/Kf)/2V{')) is related to multipath clustering in both formats, 
with V{-) denoting variance operation, respectively ll59l . 

III. Exact end-to-end Symbol Error Rate Analysis 

The end-to-end SER for the considered cooperative system can be expressed as [|48]| . Il50l 

2^-1 

Aer = E A = C.)AA = CJ (6) 

2 = 0 

where the binary vector A = [A(l), A{2), A(3),..., A{K)] of dimension (1 x K) denotes the state 
of the relay nodes in the system, with A{k) taking the binary values of 1 and 0 for successful 
and unsuccessful decoding, respectively. For the case of statistically independent channels the 
joint probability of the possible state outcomes can be represented as follows: 

K 

P(A) = P{A{1))P{A{2))P{A(3)) ■ ■■P(A{K)) = J] AA(L). (7) 

k=l 

Furthermore, = [C(l), (7(2), C(3), • • • ,C{K)] denotes different possible decoding combi¬ 
nations of the relays with e {0,2^ “ where C{k) takes the value of either 0 or 1. The 
conditional error probability P(e|A = C^) is the error probability conditioned on particular 
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decoding results at relays while P(A = C^) is the corresponding joint probability of the 
decoding outcomes. Based on the MRC method, the instantaneous SNR at the destination for 
given decoding combination, C^, can be expressed as Il50ll 

7mrc(C^) = 1^5, C{k)\aRi^^D\'^^^. (8) 

It is also recalled that the MGF for independent fading channels in DF scheme is given by fT2\ 

K 

^7MRC ) = A=f„.„(<i)nC(*^)A='7R,.„(o) (9) 

k=l 

which in the present analysis can be expressed according to |I^ eq. (6)], namely 

M (-?—\ = f __ 

"" Uin^ t)/ [ (2(h - H)i, + + H),, + ^ 7 ) 

It is noted that the above expression is particularly useful in the subsequent SER analysis. 



A. End-to-End SER for M—PSK Constellations 

1) The case of i.n.i.d rj — n fading channels: The end-to-end error probability for M—PSK 
constellations over individual rj — ^ fading link when rj, jj, and 7 in each path are not necessarily 
equal can be expressed as ITTOI eq. (5.78)] 


.Pm-psk — 


{M-l)7r 

M 


r ^ 7 ( 7 ) 

/ T9PSK 

'0 7 re e 


r-7r/2 


d7d6' = — 


M 


f gPSK \ 


TT 


nsin^^/ 


d0 + - 
vr 

—^ s 


(M-l)7r 


r /2 


I 7 


Vsin^^ 


dO 


= 2:1 


= Xii 


( 11 ) 


where f^psK = sin^(7r/M) If^ . In order to evaluate ®, we firstly need to determine the error 
probability for decoding at the destination terminal, using MRC, under given decoding outcomes 
at nodes i.e., for a given [1711 . To this end and based on the MGF approach it follows that 

/ 4:^% j^hs,D{!^{hs,D + Hs,D)ds,D + 


P(e|A = a) = - [ 

^ Jo 

K 

xHCW 


1 

TT , 
K 


k=l 

_ {M-l)n 

r M 


X 


7r/2 


- Hs,D)dS,D + 

^TRk,D^Rk,D 

‘^{hR^^D + Hr^^d)^Ir^^d + 
4/^I,Z)^5,d(2(/15,D + Hs,D)dS,D + ^f|7S,D)“^^ 


1 \ 


d9 


{2{hs,D - Hs,d)es,d + ^;t|7s,d) 


k=l 


^TRk,D^Rk,D 


-1 \ rRi,,D 


( 12 ) 


de. 






























Evidently, the derivation of an analytie solution for (fT^ is subjeet to analytie evaluation of the 
integrals Xi and Xu in elosed-form. To this end, for the ease of non-identieal fading parameters, 

f^s,D 7 ^ 7 ^ ■ ■ ■ 7 ^ Vs,d 7 ^ VRi,d 7 ^ ■ ■ ■ 7 ^ VRk,d 'ygD 7 ^ 7 ^ ''' 7 ^ lRk,D^ 

the Xi term ean be alternatively expressed as follows 

1 1 A c{k)de 


Xi = 


^ (-t , Ai NMS.d A 2 \t^S,D 11 ^ 


n 


Jhk\^^Rk•o 


where 


and 


7s,z)5'psk 


= _ 

‘^{hs,D {t} Hs,D)fJ'S,D 


(13) 


(14) 


fBik -1 _ 


7flfc,z)5'PSK 


By also setting u = sin^(0) and earrying out tedious but basie algebraie manipulations yields 

/97mrc('?PSk) /•' (1 - ^ C{k) du 


(15) 


Xi = 


27r 


1 + t 


\ MS,D 


1 +t 


MS,D 


n 


-ll+* 


MRfc.D 


1 + ^ 

t>2k 


where 

/^ 7 mrc (S'PSk) 




7|,d5'psk 


k=l 




7|;fc,_D5'psK 


tkRk,D 

(16) 

• (17) 


Importantly, eq. (fT^ ean be expressed in elosed-form in terms of [[Til eq. (7.2.4.57)], yielding 

/^7MRc(5'PSK)r(2/i5_£, + 2 Y.k=i C{k)^iR^^D + |) 


Xi = 


2v7rr(2/r5,£) + 2 C{k)fiR^^D + 1) 


(18) 


xF^ I 2/i5,D + 2 ''^C{k)iXRi^ ,D + /^5,D, /^S,D, • • • , k‘RK>Di f^Ri,D, ' ' ' , I^RkF 


k=l 

K 


; 2/is,D + 2 ^ C{k)fiRi^^D + 1; “ 1“’ 


Ai' A2' Bii ’ ’ Bik ’ i? 2 i ’ 




fc=l --- -- 

where denotes the generalized Laurieella hypergeometrie funetion of n variables fTSll . 

In the same eontext, for the Xu integral we set u = cos^(6*)/cos^(7r/M) in (O yielding 

^7MRc(^PSK)cos(7r/M) 


Xii = 


27r 


(19) 


u 


X 




1 _ cos2(7 r/M) 

1 -HAi “ 


MS.D 


1 - 


s2(7r/M) 


I+A 2 


U 


Hew- 


fc=i 


. _ cos 2 (VM) 


■dw. 



























Evidently, the above integral ean be also expressed in elosed-form in terms of the (•) funetion; 
therefore, by performing the neeessary ehange of variables and substituting in (fT9l) . one obtains 

(i?psk) 


Xii = 


7mrc ' 
TT 


K 


( 20 ) 


X I lls,D, f^Ri,D, ■ ■ ■ , f^RK,D, f^Ri,D, ' ' ' , f^RK,D, “ “ ‘^f^S,D + 2 ^ C(k) 


k=l 


cos^(7r/M) cos^(7r/M) cos^(7r/M) cos^(7r/M) cos^(7r/M) cos^(7r/M) 


1 + 1 + A2 1 + -Bii 1 + 1 + -B 21 1 + k?2K 

It is noted here that the 

Nevertheless, despite its importanee it is not unfortunately ineluded as built-in funetion in popular 
software paekages sueh as MATLAB, MATHEMATICA and MAPEE. Based on this, a simple 
MATEAB algorithm for eomputing this funetion straightforwardly is proposed in Appendix 1. 

2) The case ofi.i.d rj—fi fading channels: Eor the speeial ease of identieal relay to destination 
paths i.e., = ■ • ■ = = dR,D, VRi,d = ■ ■ ■ = r]RK,D = Vr,d, 7ri,d = ■ ■ ■ = Irj^^d = 

and thus, Bu = ■ ■ ■ = Bik = Bi, and B 21 = ■ ■ ■ = B 2 K = B 2 , equations (fTSl) and (l20l) 
ean be readily simplified to the following representations 


, cos^(7r/M) 


^ vv funetion has been studied extensively over the past deeades. 




/^7mrc (5'PSK)r ( 2/i5 D + 2^r^d Y^k=l + 


K 


2i/7rr \ 2 ^ s,d + 2 ^ r^d Yl!k=i^i^) + ^ 

dS,D-, dS,D, KfiR^D, K^r^d; 2fis,D + 2fiR^D C(k) + 1; 


A"’ 2,^s,D + 2f,„,Dy^C(A^) + -; 


k=l 


K 


k=l 


1 

1 

1 

1 \ 

211 ’ 

212 ’ 

Bf 

B2 / 


( 21 ) 


and 


M 


7MRC (5'PSk) p(5) 
^ D 


K 


TT 


/^S,D, dS,D, KfiR^D, KflR^D, - — 2 ^s,D — 2 flR^D C*(fe) 


k=l 


3 cos^(7r/M) cos^(7r/M) cos^(7r/M) cos^(7r/M) 


( 22 ) 


, cos^(7r/M) 


2 1 + Ai 1 + A2 1 + Bi 1 + B2 

respeetively. As a result, with the aid of the derived expressions for Xi and Xu, the eorresponding 
error probability for M—PSK eonstellations ean be determined by P(e|A = C^) = Xi +X 11 . 

It is reealled that the derivation of the overall SER also requires the determination of the 
deeoding probability of the relay nodes P(A = Cz). This is in faet a direet produet of the 
terms P{fs,Rk) = P{A{k) = C{k) = 0) i.e. deeoding error at the relays and (1 — P{fs,Rk)) = 


CO I CN 
















10 


P{A{k) = C{k) = 1) i.e. successful decoding at the relays, whieh, as already mentioned, is a 
pre-requisite for forwarding re-eneoded signals to the destination. Importantly, this ean be also 
determined in elosed-form with the aid of the eommonly used MGF approaeh, namely 


-7r/2 


Piris,Rk) = - 

TT , 




5'psk 


Vsin2 0 


d0 + - 

TT 

-' s_ 


(M-l)7r 

M 




'tt /2 


9psk 


Vsin^^ 


de 


(23) 


= ^X22 

whieh with the aid of (fTOl) ean be equivalently re-written as follows 


P{A{k) = C{k) = 0) = - 


TT 


(‘^ Ps,Rhs,Ru {‘^{hs,Rk + Hs,R^)^is,Rk + 
- Hs,Rk)l^S,Ru + 


-1 \ 


de 


1 

TT 


(M-l)7r 


" [‘^Ps,R^hs,Rk -f Hs,RjfiS,Rk + 


G/2 


-1 \ 


(2(/is,Rfc - Hs,R,,)9.s,Rk + ^^7s.rJ ' 

Notably, the integrals in (l24l) have similar algebraie form to Xi and Xu sinee the differenee is 
the absenee of terms. As a result, the following elosed-form expressions ean be dedueed 

/^7S,H;,(5'PSK)r(2/i5^i^j^ -f 2) ^(2) 


Xo = 


2 ^/n^ {2 fis,Rk + 1) 


-Pd^ ( + 2? k'S,R^,, k'S,R^,; ‘^k^S,Rk + 1; 


(25) 


and 


X 22 — 


where 


and 


(5'psk) 


TT 


^(3) 


1 1 3 cos^(7r/M) cos^(7r/M) ^ 


2I k'SyRkI k'SyRkI 2 ‘^k'S,Rk 1 2' \ (J 


1 + C*2 


, cos^(7r/M) 


{S} = 


7s,iJ;,5'PSK 




/^s,Rk{g\ 


pskJ = 


4/^1, ^S,Rk ) ^ 


',Rk 
'y‘s,Rk9pSK 


(26) 

(27) 

(28) 


Therefore, the P^r for M—PSK is dedueed by inserting P(e|A = C 2 ) and P(A = C^) in dH). 
To the best of the authors’ knowledge, the derived analytie expressions are novel. 


B. End-to-End SER for M—QAM Constellations 

Having derived the SER over i.n.i.d and i.i.d 9 — jj, fading ehannels for the ease of M—PSK 
modulations, this subseetion is devoted to the derivation of the eorresponding SER for the ease 
of M—QAM eonstellations. 
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1) The case of i.n.i.d rj—jj fading channels: In the case of independent but not necessarily 
identically distributed rj — ^ fading channels and based on fTOl eq. (9.21)], it follows that 


Pm- 


M-QAM 


AC 


71 


f7r/2 


M 


7mrc 




d'K j A. 


Vsin 


2^; 


TT 




^sin 


^eJ 


(29) 


= I4 


where (?qam = 3/2(M — 1) and C = {1 — I/s/M) irTTll . Therefore, the probability of decoding 
error using MRC under given decoding results at the relays can be expressed as follows: 


P( I A _ p \ f 4/^I,d^5,d(2(/is,z) + Hs,D)dS,D + f^75,Z)) ^ 

“ vr io V ms,n - Hs,D)ds,D + i^ 7 s,z)) 


-1\ MS,D 


(30) 


X 


J-J I '^TR^:,DhR,^,D + Hr^^d)pr^^d + ^ 


k=l 




HR^,D)dRk,D + ^^l'Rk,D) 


AC^ /•’r/4 / 4^2 ^^(2(/i5^^ + Hs,D)ds,D + 77;^7s,d) 


TT 

K 


tJ'S,D 


>0 \ mhs,D - Hs,D)ds,D + f^^7S,Z)) 


Y^C{k) f + Hr^^d)pr^^d + f§^7R^^,D) ^ 


X /u 

fll V (2(^Rfc,D - HR^^D)dRk,D + f^7i?fe,D) 

Evidently, the derivation of a closed-form expression for (l30l) is subject to analytic evaluation 
of the above two integrals that correspond to X3 and X4. Following the same methodology as in 
the case of M—PSK constellations and using the change of variable of u = sin^ 0, one obtains 

/97mhc('?QAm) /■Xl-M)-iu2(MS,z,+ELiC(A:);.«„,)-l _ 


X 3 = 


, k-S,D 
1 + ^ (1 + 


A 2 J 


MS,D 


n 


C(k) du 


=1(1 + 


Blk ) 




^■Rk,D 


(31) 


where the parameters Ai,A 2 ,Bik and i? 2 fc are determined by substituting gpsK with (?qam- To 
this effect and after some algebraic manipulations, the following analytic expression is deduced 

/^7MRC (5'QAM)\/^r (2fIs,D + 2 Y,k=l C{k)fiR^^D + I 


X 3 = 


2r (2fis,D + 2 C{k)fiR^^D + 1) 


K 


(32) 


X ( 2/rs,D + 2 C{k)pR^, ,D + 77] dS,D, dS,D-, dRi,D, ' ' ' , dRK,D, dRi,D-, ' ' ' , TRk,D 


k=l 

K 


1 1 

; 2/i5,D + 2 ^ C{k)pRf^^D + 1; 

k=i ^ 2 


B 


11 


B 


IK 


B^ 


21 


B. 


2,K 
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Likewise, the X 4 integral has the same integrand but a different upper limit of integration. Thus, by 
following the same methodology and setting y = 2u, yields the following elosed-form expression 

/^7mrc (5'QAM)r(2/i5,D + 2 Ylk=i + |) 


X 4 = 


r(2/i5,D + 2 Ef=i C{k)yiR,^D + Cik),a„o)+l 


K 


X ( 2 / 25 ,D + 2 C{k)fiRf, ,D + /^5,D5 • • • 5 IJ‘Rk,Dj , f^RK,D] 


k=l 


K 


; 2/i5,D + 2 ^ C{k)fiR^^D + 

where 


2 

1 1 


(33) 

1 

2 


fc=i 


/^7MRC (5'qam) — 


2 ’ 2 Ai’ 2^2 ’ 2X11’ ’ 2 Bik' 2B21 

/ a,, 2 t. \^^S,D K /a.. 2 <1.2 tt2 \\ 

/ 4/i5,D^S,Z) \ 


2 X 2^’2 


y ' y /i;=i 




_ 2 2 

Tk/cDS'qam 


(34) 


Notably, the X 3 and X 4 terms are also expressed in terms of the generalized Laurieella function. 

2) The case of i.i.d rj—fi fading channels: In this simplified scenario, the corresponding 
solutions for X 3 and X 4 for the case of M—QAM constellations can be derived by following the 
same methodology as in the case of M—PSK modulations. To this end, after a necessary change 
of variables and long but basic algebraic manipulations, it immediately follows that 


X3i,i,d 


Amrc (5'QAM)\/7rr i2y,s,D + 2/25,_d C'(fc) + I 


( 2 /i 5 ,D + 2 fiR^D ^ c{k) + 


K 


2 r (^2fis,D + 2 / 2 _r,d C{k) + 1 

1 

ds,D, ds,D, KdR,D, KfiR^D] 2fxs,D + 2 


k=l 


K 


k=l 


and 


Amrc (gQAM)r(2/i5'D + c{k) + 

r(2f.s,D + 2 ah,d Efci C(A^) + |)22 *‘ao+2m,„ Ef., cml 

1 ^ 3 

MS,D, dS,D, KfiR^R, Kf^R^R, 2/25,5) + 2/15;^,,5 ) EcW + 7 

k=l 


Likewise, the corresponding decoding error probability at the relay nodes can be expressed as 


3 

+ 2 ’“ 

1 1 

1 1 \ 

A’ A’ 

1 ? 

1 

K 

1 



(35) 

■C 

K 

2fis,D + 2fiR^j:, ''^^C{k) - 
, k=l 

1 

1 

1 1 1 

’ 2 Ai 

’ 2 A 2 ’ 

2 X 1 ’ 2 X 2 ’ 2 



(36) 

= a) 

= ACI^/tt 

-AC'^IJtt. 


1 

2 


Phs,Rk) = 


AC 


M 


( fi'QAM \ 


TT 


Vsin^^y 


de 


4(^2 r-K/i 


M 


( fi'QAM \ 


TT 


Vsin^^y 


de. 


(37) 
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Based on the approaeh in (1241) and given the non-arbitrary limits of integration, the following 
exaet elosed-form expressions are dedueed for the integrals X 5 and Ig: 


/^7S,flfe(5'QAM)\/^r(2/r5 Rj^ Y,k=l + \) 


X F, 


( 2 ) 


D 


11 

^ C{k) + k'S,R^; ‘^k‘S,Rk + 1; 

k=l ^ ^ 


(38) 


and 

Xg 


(fi'QAM)r(2/i5^R^ C{k) + i) 
r(2/X5,rt, Ef=i C{k) + 


X F; 


(3) 




1 1 3 3 

2/iS,Rfc ^ C(/c) + /i5,Rfe, /^5,Rfe, 2/iS,Rfe Ec(a,) + »; 


A:=l 


fe=l 


2 ’ 2CF 2CF2 


respeetively, where 


/^7S.flfe (5 'qam) - 


^Fs,R^ihg,Rk ^S,Rk) 

_ 2 2 

^S.itfcS'QAM 




(39) 


(40) 


whereas and C 2 are obtained by substituting ^fpsK with (?qam in (1271) . Therefore, using 
and (l39l) . the eorresponding deeoding error probability ean be readily expressed as 

X( 7 sa) = P{A{k) = C{k) = 0) = —X 5 - ^Xg. (41) 

71 71 

Thus, a elosed-form expression for the SER is dedueed by substituting P(A = C^) in along 
with the eorresponding P(e|A = C 2 ) and ean be eomputed using the algorithm in Appendix I. 


C. Simple Asymptotic Expressions 

The derivation of asymptotie expressions typieally leads to useful insights on the impaet of 
the involved parameters on the system performanee. This is also the ease in the present analysis 
as simple elosed-form asymptotie expressions are derived for high SNR values. To this end, the 
MGF of ?7 — /i distribution ean be aeeurately approximated as 

4/i^h \ ^ /4/i^h^ ^ 


M 


7»7-M 


9 


{2{h - H)p + -^7)(2(h + H)p + 


g 2 j 2 


sin^'^(0). 


(42) 


Based on this, the eonditional error probability P(e|A = C^) ean be approximated as follows: 
P(e|A = C.) « ( ) p[ a„.,c(C7 (1 


9 ls,D 


k=l 


9 lRk,D 


(43) 
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where Ar^^]:){Cz) for M—PSK constellations is given by 


(M-l)7r 


Ar^,d{C^) - - 

71 




(44) 


which can be equivalently re-written as follows: 
1 




TT 


/2 1 


(M-l)7r 

M 


TT 


t/2 




I? 


Xs 


(45) 


The derivation of a closed-form expression for (1431) is subject to analytic evaluation of the above 
trigonometric integrals. Hence, setting u = sin^( 6 *) and a = 2{fj.s,D + 'l2k=i^i^)f''Rk,D) yields 


X7= :r 


U 2 


-dw 


(46) 


(1 -m) 

which can be expressed in closed-form according to |[69l eq. (3.191.1)]. To this end, by perform 
ing the necessary variable transformation and after basic algebraic manipulations one obtains 

^/%T{2l^s,D + 2 Ylk=i C{k)fj.R^,D + 1) 


l7 = 


2'T{2^s,d + 2 Ylk=i Cik)^R^^D + 1 ) 

Likewise, by setting u = cos^{9)/ cos^( 7 r/M) to the second integral in (l44l) . it follows that 

cos(7r/M) 


(47) 


Xr = 


U 2 


/o 


(l — M cos' 


iil))" 


-du 


(48) 


which can be expressed in closed-form in terms of the Gauss hypergeometric function, yielding 

K 


Xg - cos ( — ) 2 X 1 j - - 2^s,d + 2^C'(/c)/iRj^,£), 


1 3 


MJ 


k=l 


2 '2 '”='(s) 


(49) 


To this effect, Ari^^d{Cz) can be expressed in closed-form by Ari^^d{Cz) = -t-Xg/vr. 
In the same context, for the case of M—QAM constellations one obtains 


Ar^,d{Cz) - 


AC 


sm 


TT 


(5l)]4Gs,D+Ef=iC'(fe)A‘%,D)^5) — 


4(^2 /•7r/4 


sm 


TT 


(^Q'^YCs,D+Y.k=l C{k)lkR^,D)f^Q _ 


Xg 


— 

Xio 


(50) 


It is evident that the integrals X 7 and Xg have the same algebraic forms. Hence, it follows that 

\/^r(2(/i5,£) -f X]f=i c{k)ij.R^^D) +1) 


Xg = 


‘2T{2{fis,D + Ylk=i C{k)^R^^D) + 1 ) 


(51) 
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for the first integral, while and upon setting u = 2 sin^( 0 ) in the seeond integral, one obtains 

2-(2ms,o+2) 2F1 2 fis,D + 2 Ylk=i C { k ) nR^^,D + 2/is,D + 2 Ylk=i C { k ) fiR^^,D + |; 


2^10 — 


i'll-! 


2‘^T.k-iC{k)mk,D Y(2fis,D + 2 J2k=i C{k)fiRi^^D + |)[r(2/i5^B + 2 X]fc=i C{k)^R^^D + 2)] 

(52) 


Therefore, with the aid of (|5T]) and (l52l) . a elosed-form expression for Ar^^r){Cz) for the ease 
of square M—QAM eonstellations is given by Ar^^j:){C;^) = ACX^/ti — dC^Xio/vr. 

The deeoding error probability of the relays also in the high SNR regime ean be expressed as 


P{A{k) = C{k) = 0) ^ As,r, 


'^k^‘s,Rk^S,Rk 

9^%Rk 


k'S.Rf, 


and 


P{A{k) = C{k) = 1) ^ 1 - As,r, 


'^k'S,Rk^S,Rk'\^^’^'" 


a^^lRk 

where Ar^r^ for M—PSK and M—QAM eonstellations is given by 

1 


^S,Rk - 


(M-l)7r 

M 


TT 




and 


AC AC"^ rA 

As^r^ = — / sin^'^S'^ {e)de -/ sm 

^ Jo ^ Jo 


4ms,H t. 


{e)de 


(53) 


(54) 


(55) 


(56) 


respeetively. Exaet elosed-form expressions for As^Rf. for both modulation sehemes ean be 
obtained by following the same methodology and proeedure as in the previous ease yielding 


^ _ r{2fis,R, + l) ^ cos(7r/M) ^ fl ^ 13 

20Fr(2/r5,R, + 1 ) ^ tt ^"[2 ^^ 3 ,Rk, cos {it/ ) 

for the ease of M—PSK and 

2Cr{2fis,Rk + I) C^'^{2fis,Rk + |) 2 X 1 (4, 2fis,Rk + i; 2/i5,ijfe + |; I) 


Ac — 


(57) 


(58) 


\/^r(2/i5,R^ + 1) T{2fis,Rk + 1) 

for the ease of M—QAM. It is noted here that at suffieiently high SNR, the probability P{A{k) = 
C{k) = 1) is elearly smaller than unity and thus, 1 - As,R^{{A^l^p.^^hs,R ^)/~ 1. 
To this effeet and by substituting (l43l) and (15^ into dH), the SER of the multi-relay regenerative 
system over generalized fading ehannels at the high SNR regime ean be expressed as follows: 


pD ^ 
-^SER ~ 


2 ^YYAr^A^A[ 


9As,d 


z=0 k=l 


k=l 


9As,Rk 


( 59 ) 


where g eorresponds to ^fpsK = sin^( 7 r/M) or gqAu = 3/(2(M — 1 )) aeeording to (fTTl) and 
respeetively, depending on the seleeted modulation seheme. 
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D. Amount of Fading 


It is recalled that the amount of fading (AoF) is a useful metric for quantifying the fading 
severity in the communication scenarios and is defined according to fTOl eq. (1.27)], namely 

Var(7MRc) E(7 mrc) “ (E(7mrc))^ 


AoF = 


(60) 


(E(7mrc))^ (E(7mrc))^ 

The moment of the 7mrc under the decode-and-forward strategy can be determined by II72II 

K 


Fn = (- 1 )' 






k=l 


(61) 


J s=0 


Based on this, the first two moments in (IMl) are obtained for n = 1 and n = 2, namely, 
E(7 mrc) = 5M^^^(s)/(9s|«=o and E(7 ^rc) = '9^M.^mrc(s)/'9^s|^=o- To this effect and recalling 
(l 8 ])—(fTOl) as well as carrying out long but basic algebraic manipulations, the corresponding 
minimum AoF, if all relays decode correctly, can be represented as follows: 


AoF 


Fs,d{fs,d + ~ 2(52/ig,j) — 25^KtLR^D 


Fr,d{Kfr,d + l ) d 4 + ‘^Fs,DFR,DSi5i 


- 1 
(62) 


where — 7 s,d^s,d/ Fs,D{hs j^ Efl^,), S 2 — i‘s,dI^F s,Di^‘s,D Hs,d)i ^3 — i\dI‘^FR,D i^‘R,D 
-Hi d) aud (54 — 'yR^D^R,D/FR,D{h‘ji j^ Hln)- By recalling that the 7 — /r model includes the 
Nakagami—m, Hoyt and Rayleigh distributions, the AoF for these cases can be readily deduced. 


IV. Optimum Power Allocation 

This section is devoted to the derivation of the OPA strategy that minimizes the derived 
asymptotic SER of the considered regenerative system subject to the sum-power constraint of a 
power budget P. Since the derived SER expression in (l59l) is a function of the power allocated 
at the source and the K relays of the system, the available power should be allocated optimally 
in order to enhance the end-to-end quality of the transmission. Based on this, the corresponding 
non-linear optimization problem can be formulated as follows: 

Uopt = argminPg^j^ 
a 

^ (63) 

Subject to : Uo + ^ RRk = «o > 0, > 0 

k=l 

where a = [oq, ohi, 0^21 • • • > denotes the relative power allocation vector. Importantly, the 
cost function in (1631) is convex in the parameters oq and over the feasible set defined by 
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linear power ratio constraints. The corresponding proof is provided in Appendix II. To this effect 
and following the definitions in ITTSII . the Lagrangian of this optimization problem is given by 

K 


K 


L +AI oo+ 


aR, 




(64) 


k=l 


k=l 


where A denotes the Lagrange multiplier satisfying the power constraint whereas /iq and /i^ 
parameters serve as slack variables. The nonlinear optimization problem in (l6^ can be solved 
using e.g. a line search method. However, to develop some insights for the power allocation 
policy we apply the Karush-Kuhn-Tucker (KKT) conditions for minimizing the corresponding 
SER fTSll . To this end, it follows that all yUfc and /iq parameters are zero while the following 
derivatives form the necessary condition for maximizing the performance of the system 






{l< k< K) 


(65) 


dao dttR^ 

In order to obtain feasible relations between optimal powers of the cooperating nodes, we employ 
the asymptotic SER in (I59l) . By re-writing this in terms of the power ratios allocated at the 
transmitting nodes it follows that 




D 

SER 





^k^R,D^Rk,D^0 


2 \ 


z=0 k=l 




^f^S,Rk^S,RkNo 

O^09^^S,RkP^ 


2 \ 


( 66 ) 


In order to derive an optimal power allocation policy for the DE protocol, we initially restrict 
our scenario to A' = 1 and K = 2 relay nodes and then we generalize the result for K relays. 

1) K = 1 Scenario: The possible decoding sets in this case are Co = 0 i.e., the relay is 
unable to decode correctly, and Ci = 1 i.e. the relay decodes successfully. Thus, using (l6^ and 
neglecting the constant term outside the summation after factoring out oq, it follows that 

, K2 


mm 


2(MS,D+A‘S,fll) 

Subject to : Oo + Rr^ = 1 


(67) 


where 


Ki = As,dA 


S,DCiS,Ri 


^i^S,Ri ^S,RiNq 


2 \ 


( 68 ) 
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and 


K, = A 


Ri,D 




2 \ ^«l.D 


(69) 


Next, we apply the necessary condition in (1651) to determine the relation between optimal power 
allocations in the two nodes. Thus, the first derivative of with respect to oq is given by 

^^SER + fJ^S^Ri)^! ‘^fJ‘S,DK2 


dan 


2(MS,i3+MS,Hi)+l 


a, 


2^3,0+! 


0 


Ri 


whereas the first derivative with respect to aR^ is given by 




da 


Ri 


a, 


‘^^^Rl,DK2 

'iR-S,D 2/ii{j^£) + l‘ 


(70) 


(71) 


0 


a 


Ri 


Equating the above two relations and after some long but basic rearrangements, it follows that 

/ I \ ‘^RS,Ri / 


K. 


2 


(72) 


It is noticed that the left-hand side of (1721) depends only on the channel parameters. As a result, 
this term is always positive which implies directly a power policy of ao^R^^R > aR^^j:,^s,D- 
This relation is also in agreement with the power allocation over Rayleigh, Nakagami—m, and 
Hoyt (Nakagami—g) fading distributions under the special cases given in section III. 

2) K = 2 Scenario: In this case, the following four scenarios are valid: Cq = (0,0) when 
the two relays decode incorrectly; Ci = ( 0 , 1 ) when the first relay decodes incorrectly and the 
second relay decodes successfully; C 2 = ( 1 , 0 ) when the first relay decodes successfully and the 
second relay decodes incorrectly; and C 3 = (1,1) when both relays decode successfully. Based 
on this, the corresponding optimization problem can be expressed as follows 


mm 


76.8 


+ 


K. 


+ 


K. 


+ 


Kn 


‘i{RS,D+rS,Ri+RS,R2) ArS,D+rS,R-d 2/^fl2,D ArS,D+rS,R2) ‘^rRi,D ‘2RS,D ‘^rR2,D '^IJ'Ri.D 


"^2 


^Ri 


‘'R 2 


"Hi 


(73) 


Subject to : ao+aR^+aR^ = 1 


where K^, 764 , and Kq relate to the channel parameters, which are not affecting the sign of 
the derivatives in any case. To this effect, the derivative of P^r with respect to aR-^ is given by 

2KQjJ.R^,D 


^PsER 


da 


Hi 


2(ms,d-I-MS,H2) 2aihj,_d-|-1 2/is,_d 2/1^2,o 2/iHj ^0+1 


( 74 ) 


"Hi 


"H2 


"Hi 
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whereas the derivative of -P^r with respeet to is given by 

^-Ps^R ^__ 2Kq^R2,D 

«2 Cq Uq “R2 “-Rl 

With the aid of the above two representations, a power alloeation strategy ean be proposed for 
identical channel conditions as = aR^. Likewise, applying dP^^/daR^ = daR^+i, 

{2 < k < K) li is straightforwardly shown that ur^, = aR^,+i. This power assignment indicates 
that under the total power constraint and the assumed channel conditions for the regenerative 
network, power control mechanism with the relays is not essential, but the remaining power left 
from the total power budget after the source can be simply allocated uniformly between the 
nodes. This is further elaborated and largely verified in the following section. 

V. Numerical Results 

In this Section, the offered analytic results are employed in evaluating the performance of the 
considered regenerative system for different communication scenarios. To this end, the variance 
of the noise is assumed to be normalized as iVo = 1 for all considered scenarios while transmit 



Fig. 2. SER in r; — /i fading with r — 0.5 ,77 = 1 and 
IIs.D = = OdB for 4-PSK/4-QAM 

constellation with different number of relays and EPA. 



Fig. 3. SER in 77 — /r fading with r = 0.5 ,77 = 
1 and f2s,D = ^s,Rk ~ ^Rk,D = OdB for 
4-QAM/4-PSK and 16-QAM, A = 2 and EPA. 
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powers are equally or optimally alloeated to the source and the relays. It is also noted that 
the presented results are limited to Format—1 of the rj — jj, distribution but they can be readily 
extended to scenarios that correspond to Format—2 Il59ll . 

Fig. 2 illustrates the SER performance as a function of SNR for one, two and three relays 
using equal power allocation, i.e., Pq = Pr^ = P/{P + 1 ) symmetric and balanced rj — 
fading channels for 4—PSK/4—QAM constellation. Also, the f] — ^ fading parameters are /r = 
0.5 and 77 = 1 while parameters are equal to unity. It is shown that the exact results are 
bounded tightly by the corresponding asymptotic curves from moderate to high SNR values 
while a full diversity order i.e., 2, 3 and 4, can be achieved. Table II depicts the corresponding 
diversity gains computed from the slopes for the exact and asymptotic curves along with the 
direct transmission scenario, for reference, where it is assumed that Pq = P. It is observed 
that at a target SER of 10““^ the single relay system exhibits a gain of 15dB over the direct 
transmission whereas the two and three relay systems outperform the direct scenario by about 
19.5dB and 21.5dB, respectively. In the same context. Fig. 3 illustrates the exact and asymptotic 



0 5 10 15 20 25 30 35 40 

SNR 


Fig. 4. SER vs average SNR with ^ = {0.5,1,1.5}, 
ri = (0.1, 0.9} and Q,s,d = ^s,Rk = ^Rk,D = OdB for 
4-PSK/4-QAM, TV = 2 and EPA. 



Fig. 5. SER vs average SNR with /r = {0.5,1}, 773,0 = 
0.9, tls.o = ^s,Rk ~ OdB, Q.r^^d = {0,10}dB with 
diffferent 773 ,and for 4—PSK/4—QAM, K = 2 

and EPA. 
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TABLE II 

DIVERSITY GAINS FOR ONE, TWO AND THREE RELAYS USING 4—PSK/4—QAM EOR fJ. = 0.5 AND rj= 1. 


K 

Diversity Gain (Exact) 

Diversity Gain (Asymp.) 

I 

1.96 

2 

2 

2.91 

3 

3 

3.85 

4 



Fig. 6. SER in rj — fi fading with fi = {0.5,1}, rjs.D ~ 0.9 and Qs,d = ^s,r^. = OdB, = (0, lOjdB with diffferent 

Vs,Rk VRi,d for 4—PSK/4—QAM signals, K = 2 and EPA. 


results for 4—QAM/4—PSK and 16—QAM constellations for the case of two relays with equal 
power allocation over symmetric rj — jj, fading channels with /i = 0.5 and 77 = 1 as well as 
balanced links i.e., VLs^d = ^s,Rk = ^Rk,D = ^ = OdB. It is shown that the asymptotic curves 
are almost identical to the exact ones for SERs lower than around 10“^. Therefore, it becomes 
evident that in practical system designs of DF relaying at the high SNR regime, the offered 
asymptotic expressions can provide useful insights on the system performance. 

Fig. 4 illustrates the cooperation performance of 4—QAM/QPSK system in a two relay scenario 
for yu = {0.5,1,1.5} and identical channel variance of = ^s,Rk = ^Rk,D = OdB with equal 
power allocation. It is also recalled that the case of // = 0.5 corresponds to the Nakagami—g 
(Hoyt) distribution with = 7 ]. By varying the value of r], we observe the effect of the scattered- 
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Fig. 7. SER in 77 — /i fading with p = {0.5, 1}, rj = 0.5, and Q,s,d = OdB for different fls,Rk and 4—QAM/4—PSK, 

a: = 2 and EPA. 


wave power ratio on the average SER of the eonsidered regenerative system. This verifies that 
the SER is inversely proportional to rj sinee for an indieative SER of 10“^, an average gain of 
2dB is observed when r] inereases from 0.1 to 0.9 for all values of /i. 

Eurthermore, average gains of 4dB and 1.75dB are obtained as /r inereases from 0.5 to 1 
and from 1 to 1.5, respeetively. Eikewise, Eig. 5, demonstrates the SER performanee in i.n.i.d 
?7 — /i fading ehannels for 4—QAM/4—PSK constellations for the case of two relays with equal 
power allocation. It is assumed that rjs^D = 0.9, ^s,d = ^s,Rk = OdB, = {0, lOjdB, 

/ 45 ,d = l^s,Rk = l^Rt,,D = = {0.5,1} with different values of 'qs,Ry, and set as 

{fls,Rk^lRi,D,'nR 2 ,D} = 10.1,0.1,0.9} and {0.9, 0.8, 0.9}, respectively. It is observed that the 
performance of the system improves substantially as 775 ,and increase as at a SER of 
10“^ almost 1.25dB and 1.75dB gains are achieved when {'r]s,Rk-:lRi,D-,VR 2 ,D} changes from 
{0.1, 0.1, 0.9} to {0.9, 0.8, 0.9} for ^r^^d = {0,10}dB, for the considered values of 

Eig. 6 illustrates the SER performance in 77 —/i fading conditions with /ir r = 0.5 and 77 = 0.1 
for 4—QAM/4—PSK constellations and balanced links of relative channel variance OdB for the 
case of two relays using equal power allocation with different jJ-s^Rk non-identical values of 
f^Rk,D- The figure shows that increasing at least one of value at a fixed or increasing 

both Hsfik /iRj,,D simultaneously can improve the cooperation performance. Indicatively, at 
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Fig. 8. SER with OPA in r] — jj, fading with rj = 0.5, Qs,d — ^s,Rk = OdB and ^r^,d = lOdB for different fj. and 
4-QAM/4-PSK and K = 2. 

SER of 10“'^ nearly 1.25dB and 1.75dB gains are observed when {^s,Rki f^R 2 } changes 
from {0.5, 0.5, 0.5} to {0.5, 0.5,1} and from {0.5, 0.5,1} to {0.5,1,1.5}, respeetively. Also, 
nearly 0.75dB and IdB gains are aehieved when {/iHi,/fR 2 } varies from {1,1} to {1,1.5} and 
then to {1.5,2} when fis,Rk = whereas a gain of 4dB is shown when both fis,Rk and 
inerease at the same time, for instanee, from {0.5, 0.5,1} to {1,1,1.5}. Likewise, Fig. 7, depiets 
the SER performanee for 4—QAM/4—PSK constellations for /r = 0.5 and 1, 77 = 0.5 and 
unbalanced channel links employing two relays with equal power allocation. It is shown that 
increasing either ^s,Rk or ^r^,,d can improve the SER and that the performance for the case 
of ^s,Rk = OdB and = lOdB is better than the reverse scenario i.e. ^s,Rk = lOdB and 

^Rk,D = OdB. For example, almost constant gains of IdB and 1.75dB are achieved when ^s,Rk = 
OdB and ^r^.,d = OdB increase to ^s,Rk =10dB and ^r^,d = OdB and from ^s,Rk =0dB and 
^Rk,D = OdB to ^s,Rk =0dB and = lOdB, respectively. This verifies that in the considered 
regenerative protocol, the overall performance improves more when increasing average power 
from the relays to the destination than from the source to the relays. 

Fig. 8, demonstrates the SER performance of the proposed optimal power allocation scheme 
for different values of the fading parameter /i, and assuming stronger channel variance from 
the relay nodes to the destination node with constant scattered-wave power ratio for the case of 
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TABLE III 

OPTIMAL TRANSMIT POWER ALLOCATIONS WITH DIFFERENT /is,AND FOR fJ.s,D = 0.5, 

rj = 0.5, fls,D = ^s,Rk = ^Rk,D = OdB using 4—QAM/4—PSK SIGNALS AT SNR = 20 dB. 




Po/P 1 Prj/p 

Po/P 1 PR 1 .R 2 /P 

0.5 

0.5 

0.6270 1 0.3730 

0.4832 1 0.2584 

1 

0.6572 1 0.3428 

0.4914 1 0.2543 

1.5 

0.6871 1 0.3129 

0.5048 1 0.2476 

1 

0.5 

0.5415 1 0.4585 

0.4006 1 0.2997 

1 

0.5302 1 0.4698 

0.3842 1 0.3079 

1.5 

0.5725 1 0.4275 

0.3971 1 0.3014 

1.5 

0.5 

0.5160 1 0.4840 

0.3724 1 0.3138 

1 

0.4652 1 0.5348 

0.3443 1 0.3278 

1.5 

0.5008 1 0.4992 

0.3424 1 0.3288 


two relays and 4— QAM/4—PSK constellations. It is also assumed that ^s,d = f^s,Rk = f^Rk,D = 

Vts^D = ^s,Rk = OdB, ^R^,D = lOdB and r] = 0.5. It is shown that when /i is small, for 
example yU = 0.5, the OPA provides small gain to the cooperation system, which, however, 
increases as /i increases. Indicatively, for a SER of 10“^ the optimal system outperforms the 
equal power allocation scenario by at least 1.5dB when /i = 0.5 and by 2.5dB and 3dB when 
yU = 1 and yU = 1.5, respectively. In addition. Fig. 9 depicts the corresponding SER for the case 
that fis,D = fJ‘S,Rk = f^Rk,D = = 0.5, ^s,D = ^s,Rk = OdB and ^r^,d = {0, lOjdB. It 

is observed that when ^s,Rk = ^Rk,D = OdB, OPA does not provide significant performance 
improvement for the considered DF network. On the contrary, the SER improves as the difference 
between i^s,Rk ^r,.,d increases. For example, it is noticed that for a SER of 10“^ and fi = 0.5, 
gains of 0.5dB and 2dB are achieved by OPA over the EPA scheme when ^r^,d = {0, lOjdB, 
respectively. Similarly, gains of 0.5dB and 3dBs are obtained when jj is increased to yU = 1.5 
while it is generally noticed that OPA is typically more effective than EPA, even in the low-SNR 
regime. 

Fig. 10 demonstrates the SER performance of the derived OPA and EPA scenarios for single, 
two and three relays over symmetric rj — jj, fading scenario, i.e., with constant fi and constant 
scattered-power ratios r] and unbalanced channel variances from source-to-relay and from relay- 
to-destination. It is assumed that fis,D = fJ‘S,Rk = l^Rk,D = I,?/ = 0.5 whereas Vts,D = ^s,Rk = 
OdB and ^r^,d = lOdB. It is shown that the OPA strategy clearly outperforms its EPA counterpart 
since the gain for a SER of 10“"^ is 2dB, 2.5dB and 2.75dB for one, two and three relays, 
respectively. The characteristics of the OPA strategy are further analyzed with the aid of Tables 
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Fig. 9. SER with OPA over rj — fi fading for jj,s,D — = fJ‘Rk,D — = 0.5, fls,D = ^s,Rk = OdB with different 

flR^,D for 4-QAM/4-PSK and K = 2. 


TABLE IV 

OPTIMAL TRANSMIT POWER ALLOCATIONS WITH DIFFERENT AND fls,D = PS.fls, = = P 

FOR?; = 0.5, Qs,d = = OdB and 4— QAM/4—PSK SIGNALS AT SNR = 20 dB. 




Po/P 1 Hr, /P 

Po/P 1 Prt,r,/p 

Po/P 1 Pr 

,R-,,R^/P 


1 

0.6270 

0.3730 

0.4832 

0.2584 

0.4036 

0.1988 

0.5 

10 

0.7968 

0.2032 

0.6974 

0.1513 

0.6328 

0.1224 


100 

0.9181 

0.0819 

0.8712 

0.0644 

0.8371 

0.0543 


1 

0.5925 

0.4075 

0.4368 

0.2816 

0.3520 

0.2160 

1 

10 

0.8316 

0.1684 

0.7343 

0.1328 

0.6658 

0.1114 


100 

0.9557 

0.0443 

0.9247 

0.0376 

0.8995 

0.0335 


1 

0.5735 

0.4265 

0.4131 

0.2935 

0.3274 

0.2242 

1.5 

10 

0.8496 

0.1504 

0.7549 

0.1226 

0.6850 

0.1050 


100 

0.9683 

0.0317 

0.9439 

0.0280 

0.9232 

0.0256 


Ill, IV and V, which depict the optimal power ratios allocated to the souree and the relay-nodes in 
terms of Pq/P and ^Rk /P. In ease of multiple relays, the relays are assigned with equal powers 
(Pri^/P = Pri^+i/P). The power ratios alloeated to the souree and relay nodes for asymmetrie 
and balaneed channel conditions are tabulated indieating that the numerical values are in tight 
agreement with the formulated power strategy in section IV. Table IV also corresponds to the 
ease that the souree is assigned with high proportion of power when ^r^,d is larger than ^s,Rk^ 
i.e., unbalaneed ehannel links whereas Table V verifies that the optimal power alloeation scheme 
is dependent upon the considered modulation seheme. 
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TABLE V 

OPTIMAL TRANSMIT POWER ALLOCATIONS WITH DIFFERENT MODULATIONS AND 

MS.D = MS,fife = Mflfc ,D = fj.,v = 0.5,ns,D = i^s,Rk = ^Rk ,D = OdB FOR two relays at SNR = 20 dB. 



A-PSK 

16-P5fC 

16-QAM 

M 

Pq/P I -Pr, ,r,/P 

Pq/P I PRi ,R'?JP 

Po/P|Pr, ,R,/P 

0.5 

0.4832 1 0.2584 

0.4932 1 0.2534 

0.5113 1 0.2443 

1 

0.4368 1 0.2816 

0.4392 1 0.2804 

0.4572 1 0.2714 

1.5 

0.4130 1 0.2935 

0.4138 1 0.2931 

0.4287 1 0.2857 



Fig. 10. SER with OPA over rj — jj, fading for — 1, rj = 0.5, fls.D = ^s,Rk ~ OdB and Qr^^d = lOdB for 4—PSK/4—QAM 
and different number of relays. 


VI. Conclusion 

In this paper, we analyzed the end-to-end performance and optimum power allocation of 
regenerative cooperative systems over generalized fading channels. Novel exact and asymptotic 
closed-form expressions for the SER assuming M—PSK and M—QAM modulated signals were 
derived over independent and identically distributed as well as independent and non-identically 
distributed channels. The derived analytic expressions were then used to draw insight of the 
different fading parameters in the generalized rf — jj, fading conditions and their impact on the 
end-to-end system performance. The offered results were subsequently employed in developing 
an optimum power allocation scheme which was shown to significantly outperform conventional 
equal power allocation strategy. It was also shown that the optimum power allocation scheme is 
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practically independent of the seattered-waves power ratio parameter from souree-to-destination, 
while it is dependent upon the number of multipath elusters as well as the seleeted modulation 
seheme and that it overall provides signifieant performanee enhaneement. 


Appendix I 

A MATLAB Algorithm for Computing the Generalized Lauricella Function 
The Generalized Lauricella function is defined by the following non-infinite single integral, 


F^\ei,bi, ■ ■ ■ ,bn,c;xi, 




r(c) 


t 


a—1 


(i-t) 


c—a—1 


-dt (76) 


r{a)r{c- a) Jo “ {1 - Xity^ ■ ■ ■ {1 - 

The numerieal evaluation of the above representation was also diseussed in [17^ Appendix V] and 
can be straightforwardly evaluated with the aid of the following proposed MATLAB algorithm: 


Function FD = Lauricella(a, bl, ... , bn, c, xl, ... , xn); 
f = gamma(c).*gamma(a).*gamma(c - a); 

Q = @(t) f.*t."'(a - 1).*(1 - t)."'(c - a - 1).*... 

(1 - xl.*t) ." (-bl) ... (1 - xn.*t) ." (-bn); 

FD = quad (Q,0,1) 

Appendix II 

Proof of Convexity of the Optimization Problem 


We provide the proof for the eonvexity of the SER expression by using (I59l) . For mathematieal 
traetability, we eonsider the proof for three relay-nodes. Based on this, the proof for larger number 
of nodes scenario follows immediately. To this end, the asymptotie SER can be expressed as 

, K2 , Ko 


R 


D 

SER 


+ 


+ 


2(AiS,D+MS,fli+MS,fl2+PS,fl3) 2{lls,D+IJ'S,Ri+I^S,R2) '^RR^,D 2 {ij.s^b+RS,Ri+R'S,R^) 




^R2 


Ka 


+ 


At 


+ 


Ko 


2(AlS.D+MS.fli) 2ri12,D 2{RS,D+^J.S,R2+^>■S,R^) 2{rs,D+RS,R2) 




"Rs 

Ay 


^Ri 


"Ri 


A« 


2(AiS,D+A‘S,fl3) 2/ifl2,D 2, rs,d 


"Rs 

(77) 


^Ri 


^R2 


"Ri 


^R2 


^Rs 


where Ai, • • • , Ag are related to the ehannel parameters. Eet /i(ao, • • • , o-rJ)-, ■ ■ ■ , fs{aor " ^ ors) 


be funetions which represent each term of For example, we assign 

/i(Ro) 2{fls,D+R-S,R-i+R'S,R2+t^S,R3)'’ ■^2(®0)RR3) 2{fls,D+R-S,Ri+P-S,R2) ’’ (78) 

Oq Oq 
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The second order derivative of /i(ao) w.r.t oq is given by 

_ 4:Ki{fis,D + f^S,Ri + /^S,K2 + f^S,R3){f^S,D + f^S,Ri + /^5,i?2 + /^S.Rs + \) 

d^ao HfJ'S,D+tJ'S,Ri+fJ'S,R2+f^S,R3+^) 

% 

The Hessian matrix of f2{ao,aR3), V^/2(a0) can be determined as follows: 


H{ao, ojjg) 


•iK2{RS,D+RS,Ri +RS,R2)il^S,D+RS,Ri +RS,R2+^) 

2(MS,o+MS,Hj+MS.fl2+l) 

“O “flg 

^K2fJ.R^,D (MS,D +MS,-Ri +MS,-R2 ) 

“O “^3 


4-R'2M-R3,-D(MS,-D+MS,-Ri+MS,-R2) 

2(MS,D+MS,flj^+fiS,H2 + 7) 2nfl3_i3 + l 
“O “H3 


4^r2 Mflg, D (/iflg, D +1) 



( 79 ) 


• ( 80 ) 


The principal minors of the matrix iT(ao, 07^3),-ffii(ao, aii;3) ^ 0,^722(0050^3) > 0 and-f7iiiT22 > 
H12H21. To this effect, the symmetric Hessian matrix H{ao, rr^) is positive semi-definite (PSD). 
Since ^ > 0 and H{ao,aR^) is PSD, i.e., f2(0.0, ciR^) h 0 , by the second order test in 

fTSll both /i(ao) and f2{oQ,aR^) functions are convex. Following the same methodology, it is 
shown that the functions /s, • • • , /s are also convex. Hence, by the sum rule of convexity fTSl 
Section ( 3 . 2 . 1 )] it follows that the total function is convex w.r.t ao,aR^,aR^ and rr^. 
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